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Could a plasma in quasi-thermal equilibrium be associated to 
the ’’orphan” TeV flares ? 
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TeV 7 -ray detections in flaring states without activity in X-rays from blazars have attracted much attention 
due to the irregularity of these ’’orphan” flares. Although the synchrotron self-Compton model has been very 
successful in explaining the spectral energy distribution and spectral variability of these sources, it has not been 
able to describe these atypical flaring events. On the other hand, an electron-positron pair plasma at the base of 
the AGN jet was proposed as the mechanism of bulk acceleration of relativistic outflows. This plasma in quasi- 
themal equilibrium called Wein fireball emits radiation at MeV-peak energies serving as target of accelerated 
protons. In this work we describe the ’’orphan” TeV flares presented in blazars 1ES 1959+650 and Mrk421 
assuming geometrical considerations in the jet and evoking the interactions of Fermi-accelerated protons and 
MeV-peak target photons coming from the Wein fireball. After describing successfully these ’’orphan” TeV 
flares, we correlate the TeV 7 -ray, neutrino and UHECR fluxes through p 7 interactions and calculate the number 
of high-energy neutrinos and UHECRs expected in IceCube/AMANDA and TA experiment, respectively. In 
addition, thermal MeV neutrinos produced mainly through electron-positron annihilation at the Wein fireball 
will be able to propagate through it. By considering two- (solar, atmospheric and accelerator parameters) and 
three-neutrino mixing, we study the resonant oscillations and estimate the neutrino flavor ratios as well as the 
number of thermal neutrinos expected on Earth. 
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I. INTRODUCTION 


Flares observed in very-high-energy (VHE) 7 -rays with absence of high activity in X-rays, are very difficult to reconcile with 
the standard synchrotron self-Compton (SSC) although it has been very successful in explaining the spectral energy distribution 
(SED) of blazars [1-3]. Although most of the flaring activities occur almost simultaneously with TeV 7 -ray and X-ray fluxes, 
observations of 1ES 1959+650 [4-6] and Mrk421 [7, 8 ] have exhibited VHE 7 -ray flares without their counterparts in X-rays, 
called ’’orphan” flares. 

Leptonic and hadronic models have been developed to explain orphan flares. A leptonic model based on geometrical consid¬ 
erations about the jet has been explored to reconcile the SSC model [9] whereas hadronic models where accelerated protons 
interact with both external photons generated by electron synchrotron radiation [10] and SSC photons at the low-energy tail [11] 
have been performed to explain this anomalous behavior in these blazars. Based on these models HE neutrino emission has been 
studied by Reimer et al. [12] and Halzen and Hooper [13]. In particular for the blazar 1ES 1959+650, Halzen and Hooper (2005) 
based on proton-proton (pp) and proton-photon (P 7 ) interactions estimated the number of events expected in Antarctic Muon 
And Neutrino Detector Array (AMANDA). They found that the neutrino rates were 1.8 (10 -3 ) events for pp (p 7 ) interactions. 

It has been widely suggested that relativistic jets of active galactic nuclei (AGN) contain electron-positron pairs produced from 
accretion disks [14-16]. Also electron-positron pair plasma has been proposed as a mechanism of bulk acceleration of relativis¬ 
tic outflows. In gamma ray burst (GRB) jets, this plasma ’’fireball” formed inside the initial scale ~10‘ cm is made of photons, 
a small amount of baryons and e ± pairs in thermal equilibrium at some MeV [17]. However, in AGN jets the fireball cannot be 
formed because the characteristic size is too large (3 r g ~ 10 14 cm) in comparison with GRB jets. Some authors found that if 
the pair plasma is expected to be optically thin to absorption but thick to scattering, the ”Wein fireball” could exist, even though 
for the size and luminosity of AGNs [18-20], Afterward, simulations with protons inside this plasma were performed by Asano 
and Takahara [21, 22], 

At the initial stage of the Wein fireball, thermal neutrinos will be mainly created by electron-positron annihilation (e + +e“-+ 
Z — > Vj + Dj). By considering a small amount of baryons, neutrinos could also be generated by processes of positron cap¬ 
ture on neutrons (e + + n —> p + u e ), electron capture on protons ( e~ + p — > n + v e ) and nucleon-nucleon bremsstrahlung 
(NN —> NN + i/j + Dj) for j = e, u, r. Taking into account that the temperature of Wein fireball is relativistic [19, 20], then 
neutrinos of 1 - 5 MeV can be produced and fractions of them will be able to go through this plasma. As known, the neutrino 
properties are modified when they propagate through a thermal medium, and although neutrino cannot couple directly to the 
magnetic field, its effect can be experimented through coupling to charged particles in the medium [23]. The resonance conver¬ 
sion of neutrino from one flavor to another due to the medium effect, known as Mikheyev-Smirnov-Wolfenstein effect [24], has 
been widely studied in the GRB fireball [25-27], 

Telescope Array (TA) experiment reported the arrival of 72 ultra-high-energy cosmic rays (UHECRs) above 57 EeV with a 
statistical significance of 5.1er. These events correspond to the period from 2008 May 11 to 2013 May 4. Assuming the error 
reported by TA experiment in the reconstructed directions, some UHECRs might be associated to the position of Mrk421 [28]. 
In addition, IceCube collaboration reported the detection of 37 extraterrestrial neutrinos at 4 a level above 30 TeV [29, 30], 
although none of them located in the direction of neither 1ES 1959+650 nor Mrk421, as shown in fig. 1. 

Because TeV 7 -ray ’’orphan” flares are very difficult to reconcile with SSC model, in this work we introduce a hadronic model 
by means of P 7 interactions to explain these atypical TeV flares registered in blazars 1ES 1959+650 and Mrk421. In this model, 
we consider some geometrical assumptions of the jet and the interactions between the MeV-peak photons coming from the Wein 
fireball and relativistic protons accelerated at the emitting region. Then, we correlate the TeV 7 -ray, neutrino and UHECR fluxes 
to calculate the number of HE neutrino and UHECR events. In addition, we study the resonance oscillations of thermal MeV 
neutrinos. The paper is arranged as follows. In Section 2 we show the dynamic model of the radiation coming from Wein 
fireball and its interactions with the protons accelerated at the emitting region. In section 3 we study the emission, production 
and oscillation of neutrinos. In Section 4 we discuss the mechanisms for accelerating UHECRs and also estimate the number of 
these events expected in the TA experiment, supposing that the proton spectrum is extended up to energies greater than 57 EeV. 
In Section 5 we describe the TeV orphan flares of the blazars 1ES 1959+650 and Mrk 421 and give a discussion on our results; a 
brief summary is given in section 6 . We hereafter use primes (unprimes) to define the quantities in a comoving (observer) frame, 
natural units (c=(i=k=l) and redshifts z~ 0 . 


II. ORPHAN TEV 7-RAY EMISSION 


Different hadronic models have been considered to explain TeV 7 -ray observations presented in blazars [31-34]. In those 
models, SEDs are described in terms of co-accelerated electrons and protons at the emitting region. In this hadronic model, 
we describe the TeV 7 -ray emission through 7 r° decay products generated in the interactions of accelerated protons and seed 
photons coming from the Wein fireball, as shown in fig. 2. 
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II.1. MeV radiation from the Wein fireball 


The Wein fireball connects the base of the jet with the black hole (BH). We assume that at the initial state, it is formed by e ^ 1 
pairs with photons inside the initial scale r 0 = 2 r g = 4GM, being G the gravitational constant and M the BH mass. The initial 
temperature can be defined through microscopic processes at the base of the jet (Compton scattering, 77 pair production, etc). 
At the first state, photons inside the Wein fireball are at relativistic temperature. The internal energy starts to be converted into 
kinetic energy and the Wein fireball begins to expand. As a result of this expansion, temperature decreases and bulk Lorentz 
factor increases at the first state. The initial optical depth is [19, 20] 


^e,o To 

r w,o 


(i) 


where <jt = 6.65 x 10 25 cm 2 is the Thompson cross section, Tw,o = 1/y 1 — /3j ^ 0 is the initial Lorentz factor of the plasma 
and n f , 0 is the initial electron density which is given by 
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where Lj is the total luminosity of the jet, Lsdd = 27 Tm p r g /aT is the Eddington luminosity, (7 e , 0 ) = K^{1/0 O ) / K2()-/0o) — 9 0 
is the average Lorentz factor of electron thermal velocity, 9 a = T 0 /m e is the initial temperature normalized to electron mass 
(m e ), to p is the proton mass and A, is the modified Bessel function of integral order. For a steady and spherical flow, the 
conservation equations of energy and momentum can be written as 
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Here pT is the total energy density of pairs ( p e = 2m e n e ('y e )) and photons (p 1 = 3w e n 7 60 and Pt is the total pressure of pairs 
(P e = 2 m e n e 9) and photons (/L = rri f ,ri~0). The number density of electrons and photons in the Wein equilibrium, and the 
number conservation equations are given by [35] 
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respectively. Taking into account the momentum conservation law p + (pr ^p T) d j r r = 0, eqs. from (3) to (7) and following to 
Iwamoto and Takahara [19, 20], it is possible to write the evolution of the Lorentz factor, the temperature and the total number 
density of photons and pairs as 
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respectively. Eqs. (8), (9) and (10) represent the evolution of the Wein fireball in the optical thick regime. As the Wein fireball 
expands, the photon density, optical thickness and temperature decrease. At a certain radius, flow becomes optically thin (r = 1), 
then photon emission will be radiated away. Defining this radius as the photosphere radius (r=r p h), then from eqs. (1), (2), (5) 
and (10), the photon density, radius, temperature and Lorentz factor at the photosphere can be written as 
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dph = T ~ 1/3 e o , (13) 

and 

rw,ph = 'To 7 ' 3 Tiy.o , (14) 


respectively. Here L 1 = eLj is the luminosity radiated at the photosphere and e is a parameter 0 < e < 1. One can see that n 7 , 
tph, Oph, and Tp/j only depend on the initial conditions ( 0 o , r 0 , F w,o and Lj/LEdd )■ The numerical results exhibit a radiation 
centered around 5 MeV [19]. Hence, the output spectrum of simulated photons could be described by 


( (e 7 ) if e 7 < e 7b , 

de 7 \ (e 7 b) _/3 i+/3, *(e 7 ) _/3 '*, if e 7 > e 7b , 


(15) 


where f3h ~ 2, /?; ~ 1 [36] and the peak energy around e 7 & ~ 5 MeV. Also we can define the optical thickness to pair creation 
around the peak as [37] 
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where we have taken into account that the cross section of pair production reaches a maximum value close to the Thomson cross 
section. Additionally, it is very important to say that in the optically thick regimen, the angular distribution of MeV photons is 
almost isotropic, whereas in the optically thin regime it has a skew distribution. Specifically at a distance r r g , the outgoing 
photons are distributed in the range 0° < (f> p h < 60° [19]. 


11.2. Geometrical considerations and assumptions 

We consider a spherical emitting region with a uniform particle density and radius rj, located at a distance R from the BH 
and moving at relativistic speed with bulk Lorentz factor I' :j , as shown in fig. 2 (above) [38, 39]. Seed photons coming from 
the photosphere of the ’Wein’ fireball will interact with Fermi-accelerated electrons and protons injected in the emitting region. 
Taking into account R r p h, these photons (with an angular distribution; Iwamoto and Takahara [19]) will arrive and go 
through the emitting region following different paths with an angle ( <p p h ); longer paths around the center and shorter ones as the 
trajectories get farther away from the center. To estimate the distances of any path, we find the common points of the circle, 
(y — y' 0 ) 2 + {x — x' 0 ) 2 = r 2 , and the straight line, y — y 0 = in (x — Xq), through their intersections (points a and b ) (see fig. 
2 (below)). As shown in fig. 2 (below), we assume y' 0 = R — r p h, x' 0 = 0 = xq = yo = 0 and m = tt/ 2 — <f> p h- Solving this 
equation system, we obtain that the distance between two points (a and b) as a function of (f> p h can be written as 
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Photons coming from the ’Wein’ fireball will be able to go or not through the emitting region depending on their paths (distances) 
and the mean free path A 7ie = 1 /{&t n e ) inside of it. For instance, photospheric photons going through the emitting region will 
be absorbed in it if the paths are longer than the mean free path (d > A i e )\ otherwise, (d < A 7je ), they will be transmitted. 
Defining the optical depth as a function of these two quantities (d and A 7ie ) 
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it is possible to relate this angle with the electron particle density. Taking into account that a medium is said to be optically thick 
or opaque when r > 1 and optically thin or transparent when r < 1, we write the electron particle density for this transition 
(r = 1) as 


1 / 1 +tan cj> 2 ph 

2<tt y (R - r ph ) 2 - (1 + tan <fr h ) [(i? - r ph ) 2 - r 2 ] 


(19) 


The previous equation gives the values of electron density for which photospheric photons going through the emitting region 
with a particular angle could be transmitted or absorbed. 
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11.3. P 7 interactions 


Once emitted, the MeV-peak photons also interact with protons accelerated at the emitting region. Assuming a baryon content 
in this region [32, 40-43], accelerated protons lose their energies by electromagnetic and hadronic interactions. Electromagnetic 
interaction such as proton synchrotron radiation and inverse Compton will not be considered here, we only assume that protons 
will be cooled down by p 7 interactions. The optical depth of this process is 


r i 


( 20 ) 


where the photon density (n 7 ) is given by eq. (11) and er p7 is the cross section for p 7 interactions. The photopion process 
P7 —► Nit has a threshold photon energy e t h = m n + m^/2m p where the neutral and charged pion mass are ra,o = 135 MeV 
and m w o = 139.6 MeV, respectively [44]. The two main contributions to the total photopion cross section at low energies come 
from the resonance production and direct production. 

Resonance Production (A + ). The photopion production is given through the resonances A + (1.232), A + (1.700), A + (1.905) 
and A + (1.950) where the mass rri A + and Lorentzian width IA+ for all resonances are reported by Miicke et al. [45]. The 7 r° 
threshold is ~ 145 MeV, and only the jr/ —► A + —»• pn 0 is cinematically allowed. 

Direct Production. This channel exhibits the non resonant contribution to direct two-body channels consisting of outgoing 
charged pions. This channel includes the reaction pj -7 nir + , pj -7 A ++ 7 t _ , and p 7 -7 A7 t°. For the n + threshold, the direct 
pion channel p 7 —> nir + is dominant for an energy between 0.150 GeV and 0.25 GeV. 

Other contributions to the total photopion cross section (multi pion production and diffraction) are only important at high ener¬ 
gies. 


II. 3.1. 7 r° decay products 


Neutral pion decays into two photons, 7r° -7 77, and each photon carries less than 15% of the proton’s energy e 7 ~ 
A-o/2 E p = 0.15 E p . The 7 r° cooling time scale for this process is [46, 47] 
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where l p is the proton Lorentz factor, a p jr (e 7 ) ~ a vl is the cross section of pion production. The target photon spectrum 
dn 7 /de 7 is given by eq. (15) which is normalized through J 0 °° e 7 (dro 7 /<ie 7 )<ie 7 = U 1 ~ e 7 b n 7 /2T^ ph . The threshold for P 7 
interaction is computed through the minimum proton energy, which can be written as 
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where <fi is the angle of this interaction. The photopion production efficiency /„.0 = p 1 - can be defined through the dynamical 
(t' d ~ 7 ’j/Tj) and photopion time scales [48, 49] 
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where the break photopion energy is given by 
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peak and A e' pea k correspond to the energy and the width around the resonance of total cross section, respectively. Taking into 
account the proton spectrum as a simple power law 
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the photopion production efficiency (f n o), and f n o E p ( dN/dE) v dE p = e^o 7 [dN/de)^o^ de^ o >7 , then we can write the pho¬ 
topion spectrum as 
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with the proportionality constant given by 
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From eqs. (24) and (25), the proton luminosity L p ~ AttD 2 F p = AnD^E 2 can be written as 
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where A P]7 is obtained through the photopion spectrum (eq. 25). 


II.3.2. 7r decay products 


Muons and positrons/electrons are produced through charged pion decay products ( 7 ^ 

e' _ 


(25) 


(26) 


A« ± + fplvp 
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v p /v p + v e /v e ). Although muon’s lifetime i! + dec = r p + is very short with 77 + = 2.2 /j,s and m p = 105.7 MeV, muons 
could be rapidly accelerated for a short period of time in the presence of a magnetic field ( B') and radiate photons by synchrotron 
E p [31, 34, 40, 50]. After muons decay, positrons/electrons could radiate photons e 7 = 3 g^f E' 2 at 


the same place. Therefore, muons and positrons/electrons cool down in accordance with the cooling time scale characteristic, 
t'sy n i = 6Trmf/(aT B 12 E[) up to a maximum acceleration time scale given by t' syn max = 16 E [/(3 q e B'). Here q e is the 
elementary charge and the subindex i is for // and e. Hence, the break and maximum photon energies in the observed frame are 
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where we have taken into account the Lorentz factor ratios 7 , = rrif /m 2 . Assuming that Fermi-accelerated muons and 

electrons/positrons within a volume ~ 47 r ? ,3 /3 with energies 7 * rrii and maximum radiation powers P VyT nax,i — 7777 are 

well described by broken power laws 7 ““ for 7 j < 7,7 and 7 i,b 7 i < ' Q+1) for 7^5 < 7 * < 7 i t max, then the observed 

synchrotron spectrum can be written as [51, 52] 
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Here A) is the number of radiating muons and/or positrons/electrons. 


III. NEUTRINO EMISSION 

Although in the current model there are multiple places where neutrinos with different energies could be generated, only we 
are going to consider the thermal neutrinos created at the initial stage of the Wein fireball and the HE neutrinos generated by P 7 
interactions at the emitting region (see fig. 2 (above)). 


III.l. Thermal Neutrinos 

At the initial stage of the Wein fireball, thermal neutrinos are created by electron-positron annihilation (e + + e~ —> Z —► 
Vj + v 3 ). By considering a small amount of baryons, neutrinos could be generated by processes of positron capture on neutrons 
(e + +n —» p + i7 e ), electron capture on protons ( e~ +p —> n + u e ) and nucleon-nucleon bremsstrahlung (NN NN+Vj +Vj) 
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for j = e, u, t. Electron antineutrinos can be detected indirectly (in water Cherenkov detector) through the interactions with the 
positrons created by the inverse neutron decay processes (P e + p —> n + e + ). The expected event rate can be estimated by 

N ev = T p N N a V w Je <c P (K)(^Pj dE v , (29) 


where N A = 6.022 x 10 23 g _1 is the Avogadro’s number, p^r = 2/18gem -3 is the nucleons density in water [53], V w 
is the volume of the detector, cr/' p ~ 9 x 10~ 44 E^/MeV 2 cm 2 is the cross section [54, 55], T is the observation time 
and ( dN/dE )„ is the neutrino spectmm. Taking into account the relation between the luminosity L v and neutrino flux F v , 
= 47rD 2 F l/ (< E >) = 47 t D^.(E 2 dN/dE) u , then the number of events expected will be 
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Here we have averaged over the electron antineutrino energy. After thermal neutrinos are produced, they oscillate firstly in 
matter (due to Wein plasma) and secondly in vacuum on their path to Earth. 


III. 1.1. Neutrino effective potential 


As known, neutrino properties get modified when they propagate in a heat bath. A massless neutrino acquires an effective 
mass and undergoes an effective potential in the background. Because electron neutrino (v e ) interacts with electrons via both 
neutral and charged currents (CC), and muon/tau (r+ /v T ) neutrinos interact only via the neutral current (NC), v e experiments a 
different effective potential in comparison with and v T . This would induce a coherent effect in which maximal conversion 
of v e into ;/ /( (v T ) takes place even for a small intrinsic mixing angle [24]. On the other hand, although neutrino cannot couple 
directly to the magnetic field, its effect which is entangled with the matter can be undergone by means of coupling to charged 
particles in the medium. Recently, Fraija [27] derived the neutrino self-energy and effective potential up to order m at strong, 
moderate and weak magnetic field approximation as a function of temperature, chemical potential (p) and neutrino energy (E ;/ ) 
for moving neutrinos along the magnetic field. In this approach, we will use the neutrino effective potential in the weak field 
approximation, which is given by 
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(31) 


Taking into account the condition that the plasma has equal number of electrons and positrons ( N e — N e = 0), then the neutrino 
effective potential is reduced to 


Veff = - 


4\/2 Gf rrieEv 




K 0 (ai) + 



- r 


(32) 


where K, is once again the modified Bessel function of integral order i, Gf = V / 2(? 2 /8m^ / is the Fermi coupling constant, 
B c = m 2 /e is the critical magnetic field, m^ is the W-boson mass, «; = (( + l)p/(0 o m e ) and ui = (l + l)/0 o . 


III. 1.2. Resonant oscillations 

When neutrino oscillations occur in matter, a resonance could take place that would dramatically enhance the flavor mixing 
and lead to a maximal conversion from one neutrino flavor to another. This resonance depends on the effective potential and 
neutrino oscillation parameters. The equations that determine the neutrino evolution in matter for two and three flavors are 
related as follows [27, 56], 
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III. 1.2.1. Two-Neutrino Mixing. The evolution equation for neutrinos that propagate in the medium ( v e o- ;/ /( T ) is given 
by [57] 


V e ff — jg- cos 2ip sin 2ip 

iC sin 2 ^ 0 


(33) 


where 5m 2 is the mass difference, ip is the neutrino mixing angle and V e // is the neutrino effective potential (eqs. 31 and 32). 
The oscillation length for the neutrino is given by 
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and the conversion probability by 
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with oj = \J(Veff — (5m 2 /2E v ) cos 2ip) 2 + (5m 4 /4:E 2 ) sin 2 2ip. Taking into account the resonance condition 

5m 2 , 


then the resonance length ( l res ) can be written as 
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5m 2 sin 2ip 

The best fit values of the two neutrino mixing are: Solar Neutrinos: 5m 2 = (5.6t 3 '4)xl0~ 5 eV 2 and tan 2 ip = 0.427+Q'Q2g[58], 
Atmospheric Neutrinos: 5m 2 = (2.1+p®) x 10 -3 eV 2 and sin 2 2ip = I.O+qq 3 [59] and Accelerator Neutrinos: 5m 2 ~ 
0.5 eV 2 and sin 2 2ip ~ 0.0049 [60, 61]. 

III. 1.2.2. Three-Neutrino Mixing. In the three-flavor framework, the evolution equation of the neutrino system in the matter 
can be written as 


■ d v 

'di = ^ 


(38) 


where the state vector is u = (v e , i/^, v T ) T , the effective Hamiltonian is H = U ■ Hq ■ U ^ + diag(V e ff, 0,0) with Hq = 
jfj-diag(—5m 21 , 0, 5m 2 2 ), the neutrino effective potential V e ff is defined by eqs. (31) and (32) and U is the three neutrino 


mixing matrix [62-65]. The oscillation length for the neutrino is given by 
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respectively. In addition to the resonance condition, the dynamics of this transition from one flavor to another must be determined 
by adiabatic conversion which is given by 
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sin 2^13 


dVe 


eff 


dr 


-l 


> 1 . 


( 0 o r o \ 

dVeff 

l o 2 ) 

89 


-l 


> 1. 


(42) 


Combining solar, atmospheric, reactor and accelerator parameters, the best fit values of the three neutrino mixing are, 

for sin? 3 < 0.053 : A mh = (7.41±£; 2 g) x 10~ 5 eV 2 and tan 2 ip 12 = 0.446to.“g and, for sin? 3 < 0.04 : Amis = ( 2 . 1 ^ 2 ) x 10“ 3 eV 2 
and sin 2 ip 2 s = O.SOlg^gg, [58, 66], 
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III. 1.2.3. Neutrino flavor ratio expected on Earth. The probability for a neutrino to oscillate from a flavor state a to another 
flavor state j3 in its path (from the source to Earth) is 


Sal3 


1 ^ UaiUpiUajUfii Sm^ 

j>i 


( Srriij L \ 

V ) ’ 


(43) 


where Uij are the elements of the three-neutrino mixing matrix [62-65], Using the set of oscillation parameters [58, 66], the 
mixing matrix can be written as 


/ 0.817 
U = -0.505 

V 0.280 


0.545 0.191 \ 

0.513 0.694 

-0.663 0.694/ 


(44) 


Additionally, averaging the term sin ~ 0.5 [67] for distances longer than the solar system, the neutrino flavor vector at source 
(z/ e , Vt) source and Earth (v e , vv T )Earth are related through the probability matrix given by 


*4 

v T 


Earth 


/ 0.534 0.266 0.200 \ 
0.266 0.367 0.368 
V 0.200 0.368 0.432/ 



(45) 


III.2. High-energy neutrinos 


HE neutrinos are created in p 7 interactions through n7r + channel. The charged pion decays into leptons and neutrinos, 
7r ± —>- e ± -E v^/vp + C^/z^ + v & lv e . Assuming that TeV flares can be described as n° decay products, we can estimate the 
number of neutrinos associated to these flares [68-70]. For p"/ interactions, the neutrino flux, dH v /dE v = A v is related 

with the photopion flux by [see, e.g. 13, 71, and reference therein] 


dN 


dE 


1 


— E v dE v = - l — 


dN 


de 


^ de^o ^ . 


(46) 


Assuming that the spectral indices of neutrino and photopion spectra are similar a ~ cr„ [72], taking into account that each 
neutrino carries ~ 5% of the proton energy (E v ~ 1/20 E p ) [73] and also from eq. (25), we can write the normalization factors 
of HE neutrino and photopion as 




-a+2 


TeV 


-2 


(47) 


with A Pi7 given by Eq. (25). Therefore, we could infer the number of events expected through 


N ev 


TPice Na Vi / <j vN (E v ) 
J Eth 


UN 

\dE 


V 


dE v , 


(48) 


where E t h is the threshold energy, T corresponds to the observation time of the flare [13], o v n (E„) = 6.78x10 ~ 35 {E„/TeV) 0 - 363 
cm 2 is the charged current cross section [74], p tci =0.9 g cm -3 is the density of the ice and Vi is the effective volume of detector, 
then the expected number of neutrinos inferred from this flare is 


N ev 


Tp ice Na Vi 

a - 1.363 


A v (6.78 x 10 - 35 cm 2 ) 


E u ,tii \ 
TeV ) 


P 


TeV, 


with the power index 3 = —a+1.363 and A u given by eq. (47). 


(49) 


IV. ULTRA-HIGH-ENERGY COSMIC RAYS 


It has been suggested that TeV 7 -ray observations from low-redshift sources could be good candidates for studying UHECRs [75], Also 
special features in these 7 -ray observations coming from blazars favor acceleration of UHECRs in blazars [76]. In the current model we 
consider that the proton spectrum is extended up to ~ 10 20 eV energies and based on this assumption, we calculate the number of events 
expected in TA experiment. 
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IV.l. Hillas Condition 


By considering that the BH has the power to accelerate particles up to UHEs by means of Fermi processes, protons accelerated in the 
emitting region are limited by the Hillas condition [77]. Although this requirement is a necessary condition and acceleration of UHECRs in 
AGN jets [75, 76, 78], it is far from trivial (see e.g., Lemoine & Waxman 2009 for a more detailed energetics limit [79]). The Hillas criterion 
says that the maximum proton energy achieved is 

E p ,,nax =eB' r, Tj , (50) 

where B 1 is the strength of the magnetic field. Alternatively, during flaring intervals for which the apparent luminosity can achieve L* « 10 4 ' 
erg s _1 and from the equipartition magnetic field es, the maximum energy of UHECRs can be derived and written as [43, 80] 


E 


max 


1.0 x 10 21 


e \J£b L*/10 47 ergs- 1 


where $ ~ 1 is the acceleration efficiency factor. 


(51) 


IV.2. Deflections 


The magnetic fields in the Universe play important roles because UHECRs are deflected by them. UHECRs traveling from source to Earth 
are randomly deviated by galactic (Be) and extragalactic (B eg) magnetic fields. By considering a quasi-constant and homogeneous magnetic 
fields, the deflection angle due to the Bg is 


and due to B_eg can be written as [81] 


i>a tx 3.8° 


Ep : th 
57 EeV 


1 f L ° \— x 

Jo kpc 4 M G 


IpEG tX 4° 


/ E Plth y 1 

V 57 EeV ) 


f Beg 
\TnG 


Leg \ 1,/2 / l c 
100 Mpc ) \ 1 Mpc 


1/2 


(52) 


(53) 


where Lg corresponds to the distance of our Galaxy (20 kpc) and l c is the coherence length. Due to the strength of extragalactic ( Beg — 1 
nG) and galactic ( Bg — 4 /tG) magnetic fields, UHECRs are deflected (t Peg — 4° and ipG — 3.8°) between the true direction to the source, 
and the observed arrival direction, respectively. Evaluation of these deflection angles links the transient UHECR sources with the high-energy 
neutrino and 7 -ray emission. Regarding eqs. (52) and (53), it is reasonable to associate UHECRs lying within 5° of a source. 


IV.3. Expected number of events 


TA experiment located in Millard Country, Utah, is made of a scintillator surface detector (SD) array and three fluorescence detector 
(FD) stations [82]. With an area of ~ 700 m 2 , it was designed to study UHECRs with energies above 57 EeV. To estimate the number of 
UHECRs associated to "orphan” flares, we take into account the TA exposure, which for a point source is given by E t op co(5 s )/Q, where 
E top = (5) 7 x 10 2 km 2 yr, t op is the total operational time (from 2008 May 11 and 2013 May 4), ai(<5 s ) is an exposure correction factor for 
the declination of the source [83] and S2 ~ 7r. The expected number of UHECRs above an energy E Ptt h yields 

Nuhecr = F r (TA Expos.) x N p , (54) 


where F r is the fraction of propagating cosmic rays that survives over a distance > D 2 [84] and N p is calculated from the proton spectrum 
extended up to energies higher than E Ptt h (eq. 24). The expected number can be written as 


Nuhecr = F r 


zzt 0 pUj(8s) (ct 2) 

47 t (a — l)fi d% E Pt th 


(55) 


where L p is obtained from the TeV 7 -ray observations of the flaring activities (eq. 26). 


V. APPLICATION TO 1ES 1959+650 AND MRK 421 

Following Iwamoto and Takahara [19, 20], we have considered the dynamics of a plasma made of e ± pairs and photons, and generated 
by the Wein equilibrium [35]. We have obtained the evolution equations of the bulk Lorentz factor (eq. 8 ), temperature (eq. 9) and density 
of photons and pairs (eq. 10). In particular, we have computed the radius (eq. 12), temperature (eq. 13), Lorentz factor (eq. 14) and the 
photon density (eq. 11) at the photosphere as a function of the initial conditions (r 0 , Tw,o, 0 o and Lj/LEdd)- By considering the values of 
initial radius r 0 = 2r g and Lorentz factor F w,o = (3/2) 1 / 2 , we plot (see fig. 3) the initial optical depth (left-hand figure above) and photon 
density (right-hand figure above) as a function of Lj/LEdd and L 7 , respectively, and the Lorentz factor (left-hand figure below) and radius 
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(right-hand figure below) at the photosphere as a function of initial optical depth. From the panels above we take into account the range of 
values 0.1 < Lj /LEdd <30 and 0.1 LEdd < L 7 < 30 LEdd for 0 o = 1, 2, 3 to 4, and in particular for 9 = 4 we obtain that the initial optical 
depth and photon density lie in the range of 3.8 < r 0 < 990 and 1.8 x 10 10 cm -3 < n 7 < 6.75 x 10 10 cm~ 3 , respectively. In the left-hand 
figure, we can observe that r 0 behaves as an increasing function with Lj /Lea d and a decreasing function with 0 o whereas in the right-hand 
figure, photon density is an increasing function of L 1 and a decreasing function of 9 0 . From the panels below, we consider the initial optical 
depth in the range 1 < to < 3 x 10 3 to find the values of the Lorentz factors and radius in the photosphere for the initial conditions Tw, 0 =1.0, 
1.5, 2.0, and 2.5 and r 0 =1.5 r g , 2.5 r g , 3.5r g and 4.5r g . One can see that both the Lorenz factor and radius are increasing functions of initial 
optical depth and the initial conditions (Tw,o and r 0 ). In particular for rw,o=2, the Lorentz factor lies in the range 1< F\y, P h < 30 and for 
r 0 = 3.5r g the photosphere radius at 10 14 cm < r p h < 10 15,6 cm. 

Once the MeV-peak photons are released from the photosphere, they will arrive and go through the emitting region following different paths. 
By considering that each path (see fig. 4) has a length given by eq. (17), we plot the distance of these paths as a function of angle for three 
values of R = ( 8 , 10 and 12 x 10 16 cm). In this figure, one can see that the length of the paths are shorter when angle increases, achieving a 
maximum value around 30°. For instance, taking into account R = 12 x 10 16 cm, a photon arriving with an angle <j) p h = 25° only has to 
go through a distance of d = 2 x 10 16 cm. Regarding the mean free path of photons in the emitting region, we compute the optical depth 
as a function of election particle density and the angle of the photon’s paths. Taking the condition r = 1, we obtain eq. (19) and plot the 
electron density as a function of angle <f> p h, as shown in fig 5. From this figure, we can observe two zones: the transparent and opaque zones. 
For instance, a value of electron density higher than n e > 0.12 x 10 s cm -3 for (f> p h < 10°, is optically opaque, otherwise it is transparent. 
Therefore, MeV-peak photons coming from the ’Wein’ fireball could be absorbed around the center and be transmitted with angles higher 
than 25°. These photons (out of the line of view) would interact in their wake with the Fermi-accelerated protons (through p 7 interactions) 
producing n° decay products. These TeV 7 -ray photons are resealed isotropically whereas MeV-peak photons only are directed out of the line 
of view (more than 25°). In this approach multi TeV photons can be imaged by the observer, but not the MeV-peak photons. 

We have introduced a hadronic model invoking the p 7 interactions to interpret the ’’orphan” flares at TeV energies as n° decay products and 
to estimate the neutrinos and UHECRs from these events. Using the developed Monte-Carlo event generator SOPHIA (Simulations Of Photo 
Hadronic Interactions in Astrophysics; Miicke et al. [45]), we obtain the cross sections for the resonance production (jry —> A + —^ p 7 r°) 
and the direct production (pj —> nn + ), as shown in fig. 6 . From the left-hand panel, one can see that the most important contribution to the 
cross section comes from A + (1.232) with a P1 ~ 281 fi barn for e' peak ~0.31 GeV. From the right-hand panel, we can observe that the value of 
the cross section (jry —> mv + ) at t peak ~0.31 GeV is cr p7 ~ 254 /rbarn. In both cases the width around e' peak is Ae pf , ak ~ 0.21 GeV. The 
previous values of the P 7 interaction parameters will be used to compute the TeV 7 -rays from 7 r° decay products and HE neutrinos from n + 
decay products. On the other hand, we obtain the TeV 7 -ray spectrum generated by this process (eq. 25), which depends on: the parameters 
of the proton (A p and a) and seed photon (e 7 b and n 7 ) spectrum, the size of the emission region, the P 7 cross section as well as the Lorentz 
factors of accelerated protons and target photons. Using the photopion model and the method of Chi-square (\ 2 ) minimization proposed by 
Brun and Rademakers [85], we find the values of proportionality constant (A Pi7 ) and spectral index (a) through the fitting parameters; [0] and 
[1], respectively, that best describe the photopion spectrum (eq 25) [43, 86 ] 


dN 

dt 


= [ 0 ] 


(^r i+2 


-[11+3 


(56) 


It is important to say that the photopion spectrum can undergo pair production on their way to the observer which will reduce this observed 
spectrum. In this case the proportionality proton constant (A p ) increases to reproduce the observed TeV photon flux. Based on the values of 
the dynamics of the ”Wein” fireball and those values reported after describing the multiwavelength campaigns for 1ES 1959+650 (May 2006; 
Tagliaferri and et al. [87]) and Mrk 421 (from 2009 January 19 to 2009 June 1; Abdo et al. [40]), we will analyze the flares presented in both 
blazars. We consider the target photons at break energy t 1 b= 5 MeV [18-20], the values of normalization energy eo=l TeV and Lorentz factors 
Fw= 8 , 12 and 18. 


V.l. 1ES 1959+650 

The blazar 1ES 1959+650, localized at a distance of D,=210 Mpc (z=0.047), is one of the best studied high broad-line Lacertae objects 
(HBL) [ 88 , 89]. This BL Lac object has a central BH mass of ~ 2 x 10 8 Mq [90]. It was first detected in X-rays during the Slew Survey 
made by Einstein satellite’s Imaging Proportional Counter [91]. Based on the X-ray/radio versus X-ray/optical color diagram, the source was 
classified as a BL Lac object by Schachter et al. 1993 [92]. 1ES 1959+650 was further detected in X-rays with ROSAT and BeppoSAX [93]. 
In particular, a multiwavelength campaign on 1ES 1959+650 performed by Suzaku and Swift satellites and MAGIC telescope in 2006 May 
and other historical data showed the usual double-peak structure on its SED, with the first peak at an energy of ~ 1 keV and the second one, 
at hundreds of GeV. These observations with a higher flux level in the first peak than the second one were modeled by the one-zone SSC 
[87], In 2002 May, the X-ray flux increased significantly and both the Whipple [4] and High Energy Gamma Ray Astronomy (HEGRA) [94] 
collaborations subsequently confirmed an increasing 7 -ray flux as well. An interesting aspect of the source activity in 2002 was the observation 
of a so-called ’’orphan” flare, recorded on June 4 by the Whipple Collaboration [5, 6 ] which was observed in 7 -rays without any increased 
activity in other wavelength simultaneously. Within the limited statistics, a pair of hours of observations showed the same spectral shape, the 
fit of a single power law dN/dE = No x ( E/TeV)~ a to the differential energy spectrum results in No = (7.4 ± 1.3 s tat ± 0.9 sy s) x 10 -11 
TeV -1 cm -2 s _1 and a spectral index of a = 2.83 ± 0.14 stat ± 0.08 sgs . Whereas the flaring activities recorded during the month of May 
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were well described by means of SSC emission, the '’orphan flare” registered on June 4 has been interpreted in different contexts. 

Using eq. (56), we found the best set of parameters for p 7 interaction, as shown in Table 1. In addition, we plot the SED with the fitted 
parameters (see fig. 7). 



Parameter 

Symbol 

Value 

T w, P h = 8 

Proportionality constant (10 -10 erg cm -2 s -1 ) 

[0] 

-A P)7 

1.381 ± 0.231 

Spectral index 

[1] 

a 

2.716± 0.081 

Chi-square/NDF 


x 2 /NDF 

12.84/6.0 

r W,ph —12 

Proportionality constant (10 -10 erg cm -2 s -1 ) 

[0] 


1.581 ± 0.268 

Spectral index 

[1] 

a 

2.767± 0.056 

Chi-square/NDF 


x 2 /ndf 

12.45/6.0 

Tw.pfe = 18 

Proportionality constant (10 -10 erg cm _2 s _1 ) 

[0] 


1.853 ± 0.358 

Spectral index 

[1] 

a 

2.834± 0.096 

Chi-square/NDF 


x 2 /ndf 

14.66/6.0 


Table 1. Parameters obtained after fitting the TeV flare present in 1ES 1959+650 with the 7r° spectrum. 

al. [94]. Taking 
10 16 cm and B= 
and considering 

T w,ph = T, = 18 we obtained the table 2. The number of HE neutrinos reported in this table was computed with the effective volume of 
AMANDA detector (14// — 10 _3 km 3 ) [95], therefore this number corresponds to those events that could have been expected in this detector. 


From this table, we can see that for Tj = Tvia= 18, the value of spectral index is equal to that obtained by Aharonian and et 
into account the values used by Tagliaferri and et al. [87] to describe the average SED reasonably well: Tj = 18, r y = 4.66 x 
0.25 G, the initial conditions of the Wein fireball: temperature 9 0 = 1, initial radius r 0 = 4r 9 and Lorentz factor = 1.2 


Lorentz factor 

r w,ph — fj 

18 

Initial optical depth 

To 

3375.0 

Photospheric radius (10 15 cm) 

T ph 

6.01 

Photospheric temperature 

Vph 

0.07 

Ratio of total and Eddington luminosity 

Lj/L Edd 

47.34 

Photon density (10 11 cm -3 ) 

Vl 

6.26 

Proton luminosity (10 47 erg/s) 

Lp 

1.17 

Maximum escaping proton energy (10 20 eV) 

Ep,ma.x 

2.07 

Number of UHECRs 

Nuhecr 

- 

Number of HE neutrinos (10 -5 ) 


4.39 


Table 2. Calculated quantities with our model after the fit of the TeV 7-ray flare in 1ES 1959+650. 

From the previous analysis it is important to highlight that this model required a super-Eddington luminosity Lj/LEdd = 47.34. Additionally, 
as has been pointed out muons and secondary positrons/electrons are created as charged pion decay products in p 7 interactions. We look 
into the contribution of radiation released by these charged particles. First of all as shown in fig. 8, we calculate and plot the time cooling 
scales as a function of muon (left-hand figure above) and positron/electron (right-hand figure above) energies. Comparing the time scales of 
synchrotron radiation and life-time of muons, one can see that synchrotron time scale is much longer than lifetime scale for any muon energy, 
therefore muons lose negligible energy before they decay. Additionally, comparing both figures above one can see that muons cool down 
much slower than electrons/positrons, for instance for muons and electrons/positrons with energies E y = E' e ± =10 TeV the cooling times are 
4.43 x 10 13 s and 7.95 x 10 3 s, respectively. Furthermore, we calculate the positron/electron synchrotron emission at break (199.40 eV) and 
maximum (22.09 keV) energies and then plot (figure below) the proportionality constant of the electron/positron spectrum (A syn e ±) as a 
function of number density of positrons/electrons (N e ±). In this figure, one can see that for a number density of positrons/electrons 10 cm -3 , 
the electron/positron synchrotron spectrum would be A syn e ± ~ 10 -11 erg cirD 2 s _1 . The SSC flux can be roughly obtained though 
the Compton parameter Y ~ U £ /Ub is A ssc e ± ~ YA syn e ±, where electron and magnetic energy densities are U e = m e N e and 
Ub = B 2 /8n, respectively. Hence, if we consider the number density of electrons/positrons of 10 cm -3 and magnetic field B= 0.25 G, then 
the SSC flux is A sac ^_ e ± ~ 10 -14 erg cm~ 2 s _1 . 

On the other hand, target photons from the Wein fireball could also be up-scattered by electrons accelerated at emitting region. In this case 
and taking into account once again the values of electron Lorentz factors reported by Tagliaferri and et al. [87], ( 7 e ,t = 5.7 x 10 4 and 
7 e ,max = 6.0 x 10 5 ), the energies of scattered photons can be calculated by j — 2"j 2 _, e- u b, with subindex i for break and max. Considering 
the values of target photons ~ 5 MeV, we obtained that the energies of scattered photons are not in the energy range considered during the 
flare (e™ break ~ 3.24 x 10 1B eV and t": rrla;r ~ 3.60 x 10 18 eV). 
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V.2. Mrk 421 

At a distance of 134.1 Mpc, the BL Lac object Mrk 421 (z=0.03, 96) with a central BH mass of ~ 2 x 10 8 M 0 [97] is one of the closest 
and brightest sources in the extragalactic X-ray/TeV sky. This object has been widely studied from radio to 7 -ray at GeV - TeV energy range. 
Its SED presents the double-humped shape typical of blazars, with a low energy peak at energies ~ 1 keV [98] and the second peak in the 
VHE regime. As for the other blazars, the low energy peak is interpreted as synchrotron radiation from relativistic electrons within the jet; the 
existing observations do not allow for an unambiguous discrimination between leptonic and hadronic mechanisms as responsible for the HE 
emission. For instance, Abdo et al. [40] found that both a leptonic scenario (inverse Compton scattering of low energy synchrotron photons) 
and a hadronic scenario (synchrotron radiations of accelerated protons) are able to describe the Mrk 421 SED reasonably well, implying 
comparable jet powers but very different characteristics for the blazar emission site. 

Although flaring activities at TeV 7 -ray [ 8 , 99-107] with X-ray counterpart [96, 108-113] have been usually detected, atypical TeV flares 
without X-ray counterpart have occurred [7, 114]. Firstly, at around MJD 53033.4 a TeV flare was present when the X-ray flux was low. This 
X-ray flux seems to have peaked 1.5 days before the 7 -ray flux, and secondly, the higher TeV flux during the night of MJD 54589.21 was not 
accompanied by simultaneous X-ray activity. 

Using eq. (56), we found the best set of parameters for P 7 interaction, as shown in table 3. In addition, we plot the broadband SED with the 
fitted parameters (see fig. 9). 



Param. 

Symbol 

Flare 1 

Flare 2 

r w,ph — 8 

P. constant (10 -10 erg cm - 2 s -1 ) 

[ 0 ] 

Ap, 7 

5.663 ± 0.479 

3.179 ± 0.995 

Spectral index 

[ 1 ] 

a 

3.082± 0.091 

2.786± 0.045 

Chi-square/NDF 


X 2 /NDF 

6.711/6.0 

33.31/10.0 

r W,ph — 12 

P. constant (10 -10 erg cm - 2 s -1 ) 

[ 0 ] 

A 

6.138 ± 0.535 

4.078 ± 0.132 

Spectral index 

11 ] 

a 

3.143± 0.095 

3.018± 0.053 

Chi-square/NDF 


X 2 /NDF 

5.24/6.0 

36.15/10.0 

Lw,ph = 18 

P. constant (10 -10 erg cm - 2 s -1 ) 

10 ] 

A 

■^p ,7 

8.810 ± 0.854 

5.867 ± 0.189 

Spectral index 

[ 1 ] 

OL 

3.412± 0.111 

3.272± 0.053 

Chi-square/NDF 


X 2 /NDF 

3.148/6.0 

18.96/10.0 


Table 3. Parameters obtained after fitting the TeV flare presented in Mrk 421 with the 7T° spectrum. 

Requiring the values used by Abdo et al. [40] to explain the whole spectrum: T/ = 12, rj = 3.11 x 10 16 cm and B = 50 G, the initial 
conditions of Wein fireball: temperature 9 0 = 2, initial radius r 0 = 1.5 r g and Lorentz factor T w,o = 2 and considering T w,ph = T 7 = 12 
we obtained table 4. The value of HE neutrinos reported was calculated with the effective volume of Ice Cube (14// — 1 kin 8 ) [30]. 




Hare 1 

Flare 2 

Lorentz factor 

Fw,ph — Fj 

12 

12 


Initial optical depth 

T 0 

216.0 

216.0 

Photospheric radius (10 15 cm) 

Tph 

5.31 

5.31 

Photospheric temperature 

@ph 

0.33 

0.33 

Ratio of total and Eddington luminosity 

Lj / L Edd 

15.21 

15.21 

Photon density (10 11 cm -3 ) 

r/'y 

0.85 

0.85 

Proton luminosity (10 45 erg/s) 

Lp 

2.86 

8.89 

Maximum escaping proton energy (10 21 eV) 

Ep max 

18.47 

18.47 

Number of UHECRs (10 _3 ) 

Nuhecr 

0.10 

1.42 

Number of HE neutrinos 

N ev 

0.12 

0.09 


Table 4. Calculated quantities with our model after the fit of the TeV 7 -ray flare in Mrk 421. 

As pointed out for the blazar 1ES 1959+650, a super-Eddington luminosity Lj/L em = 15.21 is necessary to describe the TeV flares for Mrk 
421. By considering again the values reported by Abdo et al. [40], we plot (see fig. 10) the time cooling scales of muon (left-hand figure above) 
and positron/electron (right-hand figure below), and the proportionality constant of the positron/electron spectrum (A syn £ ± ) (figure below). 
Comparing the time scales of muons (synchrotron radiation and life-time), we can see that both scales start to be similar at E^ = 10 16 eV, 
then in this energy range muons could radiate before they decay. From figures (9) and (10) and taking into account the break photon energies 
radiating by positron/electrons (e 7 , 6 r ~ 73.64 eV and e 7 , max ~ 13.09 keV), we can see that for the number density of positrons/electrons 
greater than N e ± > 10 2 cm -8 the contribution of the synchrotron flux to the SED could be noticeable (yF„ ~ 10 -11 erg cnD 2 s^ 1 ). The 
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SSC flux for this case would be A sscn _ e ± ~ 10 -17 erg cm -2 s _1 . 

Doing a similar analysis on 1ES 1959+650, electron accelerated in the emitting region could up-scatter photons from the Wein 
fireball. In this case, from the values of the electron Lorentz factors reported by Abdo et al. [40] (7 e m in = 7 x 10 2 and 
'Je,max = 4.0 x 10 4 ), photons in the energy range ~10 TeV < e* c <1.6 xlO 4 TeV would be expected. Eventually, these TeV 
photons will undergo pair productions with the 5 MeV photons in the Wein fireball or less energetic photons around. 


V.3. Neutrino Oscillations 

We have studied the active-active neutrino process in the Wein fireball. We have used the neutrino effective potential at the 
weak field limit obtained by Fraija [27]. It depends on temperature ( 6 ), chemical potential (//), magnetic field (B) and neutrino 
energy (E„) for N e = N e (eq. 32) and N e ^ N e (eq. 31). As shown in fig. 11, we plot the neutrino effective potential at the 
weak field limit for N e = N e . Taking into account the range of values of initial temperature 1 < 0 o < 4, one can see that the 
neutrino effective potential for /x = 0.5, 5 and 50 keV lies in the range from ~ 10 -9 eV to ~ 10 25 eV. It can be seen that the 
neutrino effective potential is negative and then, due to its negativity (V e // < 0 eV), only anti-neutrinos can oscillate resonantly. 
From the resonance conditions (eqs. 36 and 40), we have obtained the contour plots of initial temperature and chemical potential 
(fig. 12). Our analysis shows that chemical potential lies in the range 0.1 eV < /x < 50 eV for solar parameters, 0.09 eV< // < 
100 eV for atmospheric parameters, and 10 eV < /x < 3.3 x 10 3 eV for three-neutrino mixing. By considering a CP-asymmetric 
7 and e ± fireball, where the excess of electrons could come from the electrons associated with the baryons within the fireball, 
we plot the neutrino effective potential at the weak field limit (eq. 31) (above) and the magnetic field contribution (below), i.e. 
by subtracting the effective potential with B=0, as a function of temperature, as shown in fig. 13. Taking into account the range 
of values of initial temperature 1 < 9 0 < 4, one can see that the neutrino effective potential for /t = 5 keV lies in the range 
from 7.6 x 1CT 9 eV to 3 x 10 ' eV (left panel) and the magnetic field contribution is the opposite as compared to the plasma 
contribution for any value of fi (right panel). From fig. 13 (left panel) you can see that the neutrino effective potential is positive 
and then, due to its positivity (V e ff > 7.6 x 1 0 9 eV), neutrinos can oscillate resonantly. From the resonance conditions (eqs. 
36 and 40), we have obtained the contour plots of initial temperature and chemical potential (fig. 12). Our analysis shows that 
chemical potential lies in the range 0.1 eV < /x < 50 eV for solar parameters, 0.19 eV< p < 95 eV for atmospheric parameters, 
0.8xl0 3 eV < /x < 4.3 x 10 4 eV for accelerator parameters and 8 eV < // < 3.2 x 10 3 eV for three-neutrino mixing. One 
can observe that initial temperature is a decreasing function of chemical potential which gradually increases as neutrino energy 
decreases. In addition, we have plotted the resonance lengths, as shown in fig. 14. As we can see the resonance lengths lie in 
the range ~ 10 4 -10 8 cm for neutrino energy in the range 1 MeV < E„ < 10 MeV. One can observe that regardless the values 
of mixing parameters and neutrino energy, the resonance lengths are much less than ~ r g , therefore neutrinos will oscillate 
resonantly many times before leaving the Wein fireball. 

Taking into account the neutrino effective potential (eq. 31) and the vacuum oscillation effects between the source and Earth (eq. 
45), we estimate the flavor ratio expected on Earth for four neutrino energies (E „ = 1, 3, 5 and 10 MeV), as shown in table 5. It 
is important to clarify that we consider a small amount of baryons, then in addition to electron-positron annihilation we take into 
account inverse beta decay and nucleonic bremsstrahlung. For the electron-positron annihilation and nucleonic bremsstrahlung 
processes, neutrinos of all flavors are created with equal probability whereas for the inverse beta decay process just electron 
neutrinos can be generated. 


Energy 

(MeV) 

4>*e : <t>^ ■■ 

r 0 — 1.0 x 10 

14 

cm 

(f) Ue : (j> Ufl : (p„ T 
r Q — 5.0 x 10 14 cm 

<t>«e : </■•> : 
r Q = 1.0 x 10 

15 cm 

6 0 =2 

1 

1.2383 : 0.8870 

: 0.8747 

1.2387 : 0.8868 : 0.8745 

1.2383 : 0.8703 

: 0.8747 

3 

1.2389 : 0.8868 

: 0.8744 

1.2390 : 0.8867 : 0.8743 

1.2388 : 0.8868 

: 0.8744 

5 

1.2391 : 0.8866 

: 0.8743 

1.2390 : 0.8867 : 0.8743 

1.2389 : 0.8867 

: 0.8744 

10 

1.2391 : 0.8866 

: 0.8743 

1.2391 : 0.8867 : 0.8743 

1.2391 : 0.8866 

: 0.8743 

e Q = 4 

i 

1.2389 : 0.8867 

: 0.8744 

1.2389 : 0.8867 : 0.8744 

1.2391 : 0.8867 

: 0.8743 

3 

1.2391 : 0.8866 

: 0.8743 

1.2391 : 0.8866 : 0.8743 

1.2391 : 0.8864 

: 0.8743 

5 

1.2391 : 0.8866 

: 0.8743 

1.2391 : 0.8866 : 0.8743 

1.2391 : 0.8866 

: 0.8743 

10 

1.2391 : 0.8866 

: 0.8743 

1.2391 : 0.8866 : 0.8743 

1.2391 : 0.8866 

: 0.8743 


Table 5. The flavor ratio of thermal neutrinos expected on Earth. For this calculation, we use the neutrino energies E = l. 3, 5 and 10 MeV, initial radius r D = 1.0 X 10 14 , 

5.0 x 10 14 and 1.0 x 10 15 cm and initial temperatures X = 2 and 4. 

In addition to the flavor ratio of thermal neutrinos, we calculate the number of events that could be expected. From eq. (30) and 
taking into account the volume of Hyper-Kamiokande experiment ~ 10 5 m 3 [115], the distance of the BE Fac objects considered 
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D z = 100 Mpc, luminosity L u ~ 10 47 erg/s and neutrino energy E v = 3 MeV, we would expected 8.14 x 10 5 events. 


VI. SUMMARY AND CONCLUSIONS 


We have considered the target photons coming from the Wein fireball, the geometry of the emitting region and the protons 
accelerated at emitting region in order to explain through p 7 interactions the ’’orphan” TeV flares presented in 1ES1950+650 
and Mrk421. Taking into account the values reported by the semianalytical description of the Wein fireball dynamics [19, 20] 
and required for the descriptions of the SEDs of 1ES1950+650 [87] and Mrk421 [40], we have naturally interpreted the TeV 
’’orphan” flare as 7 r° decay products, as shown in figs. 7 and 9. Using the values of these parameters in our model for \'w.pk = T'j 
we have computed optical depth (t q ), photospheric radius (r 7 ,/,J, photospheric temperature (9 p h), ratio of total and Eddington 
Luminosity ( Lj /E/•>/,/), photon density (//,,), proton luminosity (L p ), number of UHECRs (Nuu ecr), number of HE neutrino 
(N ev ), as shown in tables 2 and 4. From these tables one can see that in both cases, this model which has been assumed a 
spherical symmetry ’Wein’ fireball, required a super-Eddington luminosity. However, the required total luminosity becomes 
much smaller if the outflow is collimated within a small opening angle, as observed. It is important to clarify that when the jet is 
collimated to the solid angle Clj [116-118], the dynamics will not be much different from the spherical ’Wein’ fireball and then 
the total luminosity of the jet becomes smaller by a factor ilj/An. For a typical value of solid angle O ~ 10 -2 , the required 
total luminosity can be much smaller than the Eddington Luminosity. Future simulations will be needed to solve the collimation 
problem [19, 20]. 

We have showed that protons in the emitting region could be accelerated up to 2.08 x 10 20 eV and 1.84 x 10 22 eV for 1ES 
1959+650 and Mrk 421, respectively, then assuming that proton spectrum can be extended at energies greater than 57 EeV and 
taking into account the GZK limit [119, 120] we estimated the number of UHECRs expected from Mrk421 in TA experiment. 
Taking into account that 30% of proton fraction will survive to a distance of 134 Mpc [84], UHECRs cannot be expected during 
this flare as shown in the table 4. 

In addition, correlating the 7 -ray and neutrino fluxes we have computed the number of neutrinos expected on AMANDA and 
IceCube, for 1ES 1959+650 and Mrk 421, respectively. In both cases the number of neutrinos would have been 4 x 10 -5 and 
0.09 for 1ES 1959+650 and Mrk 421, which is in agreement with the non-detection of HE neutrinos in the direction of 1ES 
1959+650 and Mrk 421 (see. fig. 1). 

We have studied the active-active neutrino oscillation processes for two- and three-neutrino mixing in the Wein fireball. We have 
studied the resonance condition for a neutral plasma (eq. 32) and CP-asymmetric plasma (eq. 31). We found that for a neutral 
plasma the chemical potential lies in the range 0.1 eV < p, < 50 eV for solar parameters, 0.09 eV< // < 100 eV for atmospheric 
parameters, and 10 eV < // < 3.3 x 10 3 eV for three-neutrino mixing and for a charged-asymmetric plasma the chemical 
potential lies in the range 0.1 eV < fj, < 50 eV for solar parameters, 0.19 eV< 7 <95 eV for atmospheric parameters, 0.8x 10 3 
eV < fi < 4.3 x 10 4 eV for accelerator parameters and 8 eV < /1 < 3.2 x 10 3 eV for three-neutrino mixing. Considering the 
neutrino oscillations in a CP-asymmetric Wein fireball and in vacuum, we have computed the flavor ratio expected on Earth, 
as shown in table 5. In addition, considering a volume of the Hyper-Kamiokande experiment [115] we have computed that the 
number of thermal neutrinos expected is of the order of 10 -5 . However, if the source is of the order of Mpc with a luminosity 
> 10 47 erg/sec, a few events could be expected on this experiment. 

We have proposed that the observations of ’’orphan” TeV flares would be directly related with the numbers of thermal and HE 
neutrinos, UHECR events and flavor thermal neutrino ratios, as calculated in tables 2, 4 and 5, respectively. Although the num¬ 
bers of these events were less than one, improvements in the sensibility of detectors for cosmic rays as well as new techniques 
for detecting neutrino oscillations will allow us in the near future to confirm our estimates. 

Recently, Murase et al. (2014) studied high-energy neutrino emission in the inner jets of radio-loud AGN [49]. Considering the 
photohadronic interactions the authors explored the effects of external photon fields (broadline and dust radiation) on neutrino 
production in blazars. They found that the diffuse neutrino intensity from radio-loud AGN is dominated by blazars with 7 -ray 
luminosity of > 10 48 ergs _1 , and the arrival directions of neutrinos in the energy range ~ 1 - 100 PeV correlate with the lu¬ 
minous blazars detected by Fermi. In our model we studied high-energy neutrino emission correlated with the observations of 
’’orphan” TeV flares presented in blazars 1ES1950+650 and Mrk421, when these atypical observations are naturally explained 
through the photohadronic interactions between the Wein fireball photon field and the protons accelerated at emitting region. 
We found that although neutrinos in the energy range of ~ 1- 100 TeV are generated together with these atypical observations, a 
7 -ray luminosity of > 10 4 ' erg/sec is necessary to expect more than one event in IceCube observatory. In both models one can 
see similar conclusions, although the unusual TeV 7 -ray and neutrino correlations have been only present in two blazars. In a 
forthcoming paper we will present a more detailed neutrino estimation from the Wein fireball model. 

It is important to say that although VHE (E > 100 GeV) 7 -rays coming from these pion decay products can be attenuated through 
pair production with low energy photons from the diffuse extragalactic photon fields in the ultraviolet (UV) to far-infrared (FIR) 
(commonly referred to as extragalactic background light; EBL), the attenuation of the TeV photon flux for 1ES1959+650 and 
Mrk 421 sources which have a small redshift (z <C 0.2) is low [121], 

Finally, correlations of TeV 7 - and X-ray fluxes in flaring activities by observatories like High Altitude Water Cherenkov 
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(HAWC) which will have a sensitivity of 5 x 10“ 13 cm _2 s _1 (~ 50 mCrab) above 2 TeV [122] and satellites like Nuclear 
Spectroscopy Telescope array (NuSTAR) with a sensitivity of 2x 10 15 erg cm -2 s -1 (lx 10 14 erg cm -2 s -1 ) for an energy range 
6-10 keV (10-30 keV)[123] and Swif-XRT with a sensitivity of 3x 10 13 erg cm -2 s -1 in the range 0.3 - 10 keV[124] will pro¬ 
vide information on the amount of hadronic acceleration in blazars as well as HE neutrino fluxes. On the other hand, taking into 
account the energy range of both the 37 events detected by IceCube [29, 30] and operation of HAWC (100 GeV - 100 TeV), and 
considering the relation between the neutrinos and photons from p 7 interactions E v ~ E 1 / 2, HAWC will provide unparalleled 
sensitivity to co-relate the TeV flares and HE neutrinos. In addition, it is significant to say that NuSTAR and Swift-XRT could 
put lower limits on the synchrotron flux from secondary electrons/positrons produced by n + decay products thus confirming our 
model, as shown in figs. 8 and 10. 

Further observations are needed to understand this atypical TeV flaring activity. 
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FIG. 1. Sky-map with the location of the blazars Mrk 421 and 1ES 1959+650 (marked with red color), 72 UHECRs collected by TA experi¬ 
ments (green points) and 35 VHE neutrinos detected by IceCube (blue triangles). The red circles around the blazars have a size of 7°. 
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FIG. 2. Outline of the jet and geometrical considerations. In the figure above is showed the evolution of the Wein fireball in three stages: Initial 
state, photosphere and emitting region and in the figure below is showed the geometrical diagram of a circle (emitting region) and a straight 
line (photon path). 
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FIG. 3. Figures above: Initial optical depth (r 0 ) as a function of Lj/LEdd (left panel) and photon density n 7 as a function of L 7 (right panel) 
are plotted. In both figures we use do for 1,2 3 and 4. Figures below: Lorentz factor (Tw, P h) (left panel) and photospheric radius ( r p h ) (right 
panel) as a function of optical depth (r 0 ). 
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FIG. 4. In the left-hand figure is shown a sketch of different MeV - photon paths through the emitting region whereas in the right-hand figure, 
we plot the distance of the MeV photons as a function of angle 4>ph ■ 



FIG. 5. Electron density as a function of angle (j> p h for r = 1. In this figure, one can see that the optically thick and thin zones. 
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FIG. 6. Total cross section of the resonance production (py 


■ p7T°) (left) and the direct production (py —> nn + ) (right). 
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FIG. 7. Above: Fit of the TeV flare with our hadronic model for T =8 (dashed line in black color), 12 (dotted line in red color) and 18 
(dot-dashed line in green color). Below: SED of the multi wavelength campaign in 2006 May plus the TeV "orphan” flare in 2002 June fitted 
with our model. 
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FIG. 8. Time scale of muon and e 1 * 1 secondary synchrotron radiation presented in blazar 1ES 1959+650. Left-hand figure above: Cooling 
time scales for muons. Continuos line in green color is synchrotron radiation scale and the dotted line in blue color is decay scale. Left-hand 
figure above: Cooling time scale for electron/positron synchrotron radiation. Figure below: Proportionality constant of synchrotron radiation 
spectrum as a function of number density of electrons/positrons. 
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FIG. 9. Left-hand figure above: Fit of the first TeV flare with our hadronic model for T = 8 (dashed line in black color), 12 (dotted line in red 
color) and 18 (dot-dashed line in green color. Right-hand figure above: Fit of the second TeV flare with our hadronic model for T = 8 (dashed 
line in black color), 12 (dotted line in red color) and 18 (dot-dashed line in green color). Below: SED of the multi wavelength campaign in 
2009 plus the TeV ’’orphan” flares during the night of MJD 54589.2 fitted with our mode for Tj = 12. 
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FIG. 10. Time scale of muon and e 1 * 1 secondary synchrotron radiation presented in blazar Mrk 421. Left-hand figure above: Cooling time 
scales for muons. Continuos line in green color is synchrotron radiation scale and the dotted line in blue color is decay scale. Left-hand 
figure above: Cooling time scale for electron/positron synchrotron radiation. Figure below: Proportionality constant of synchrotron radiation 
spectrum as a function of number density of electrons/positrons. 
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FIG. 11. Neutrino effective potential in the weak magnetic field regime (left-hand figure above) and contour lines of normalized temperature 
(9 0 ), chemical potential (/r) and neutrino energy (E„) for which the resonance condition is satisfied for N e — N e = 0. We have applied the 
neutrino effective potential with magnetic field B=1 G and used the best fit values of the two; solar (right-hand figure above), atmospheric 
(left-hand figure below), and three (right-hand figure below) neutrino mixing. 
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FIG. 12. Contour lines of normalized temperature (9 0 ), chemical potential (pi) and neutrino energy (E„) for which the resonance condition is 
satisfied. We have applied the neutrino effective potential with magnetic field B=1 G and used the best fit values of the two; solar (left-hand 
figure above), atmospheric (right-hand figure above), and accelerator (left-hand figure below), and three (left-hand figure below) neutrino 
mixing. 
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FIG. 13. Neutrino effective potential in the weak magnetic field regime (left panel) and the magnetic field contribution (I4//(B = 1G) — 
Veff (B = 0)) (right panel) as a function of temperature (9 0 ). We use a neutrino energy E v = 1 MeV and chemical potential fi = 0.5, 5 and 
50 keV. 
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FIG. 14. Resonance length as a function of neutrino energy for the best fit parameters of two- (solar, atmospheric and accelerator) and 
three-neutrino mixing. 









